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Abstract 



We present a Chern-Simons (CS) theory for fractional quantum Hall states, in which 
flux attachment is followed by an enlargement of Hilbert space to include n magneto- 
plasmon oscillators and n constraints, n being the number of electrons. The simplest 
approximation yields correlated wavefunctions for fractions v = p/(2ps + 1). For 
v = l/(2s + 1) we obtain multi-quasihole wavefunctions with correct gaussian and 
normalization factors, operators that create quaisholes, and a possible Read opera- 
tor. Our composite fermion and boson operators attach flux as well as the correlation 
hole to electrons, i.e., bind them to zeros. We perform a further canonical trans- 
formation which decouples the the particles and oscillators in the long wavelength 
limit. The n constraints now act on just the particles, restricting them to the Low- 
est Landau level (LLL). Their kinetic energy vanishes In the noninter acting limit, 
they couple to the external field through a magnetic moment e/2m, carry charge 
e* = e/(2ps + 1) and a dipole moment which alone survives when v = 1/2. When a 
Coulombic interaction is turned on, a nonzero 1/m* ~ (2^) 2 / 3 e 2 £o/6 arises. Density- 
density correlation functions and Hall conductivity, additive over the particles and 
oscillators, are calculated within an approximation for v = 1/2. The formalism 
makes precise and unifies many prevalent notions regarding composite fermions and 
bosons, quasiparticles, cyclotronic physics, m and m* . 

I. INTRODUCTION 

The spectacular discovery of the Fractional Quantum Hall Effect (FQHE) by Tsui, 
Stormer and Gossard! has sparked a remarkable amount of activity among theorists and 
experimentalists!'!. The theoretical efforts fall into two broad categories: those that seek 
trial wavefunctions that capture the essential physics, for example references!~0, and those 
that attempt to construct a theory that starts with electrons in a magnetic field subject to 
Coulomb interactions and try to work their way towards the experiments through a series 
of approximations. Chern-Simons (CS)0 field theories play a central role in the second 
category^H^. The present article, also of this genre, is an elaboration of a recent letterS. 
It relies on its predecessors, extends them, unifies many miscellaneous ideas that have been 
in the air and exposes the underlying physics in a particularly transparent way. 

We follow the standard notation in which the electrons, n in number, have a bare (band) 
mass m, charge — e, a complex coordinate z = x + iy, move in unit area in the x — y plane, 
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and are in a field B = —2B pointing down the z— axis. We set % = c = 1. The inverse 
filling fraction 

^ = ^- (1) 
2ixn 

gives the number of flux quanta per electron. In the cases of interest, v < 1, it will be 
assumed all spins are polarized. 

Let us first neglect interactions. All the electrons will populate the Lowest Landau Level 
(LLL). The ground state will be macroscopically degenerate and the first nonzero energy 
excitation will be at the cyclotron frequency 

u c = — . 2 

m 

The effective mass of the particles will be infinite: 

A = (3) 

m* 

reflecting the quenching of kinetic energy. 

Let us now turn on a Coulomb interaction, weak in the sense 

e 2 /l « uj c (4) 

where the cyclotron length 

k = (eB)- 1/2 (5) 

is a measure of the typical interparticle spacing and the size of the orbits, and where the 
dielectric constant e has been suppressed. One may reasonably expect that in the limit 
m — > all higher Landau levels and any reference to m will disappear from the picture; that 
the degenerate manifold of the LLL states will broaden out to a width controlled by the 
effective mass 1/m* ~ e 2 /o, of Coulombic origin. Thus consider Laughlin's wavefunctionsB 
at v = 1/3 for the ground state, and state with a quasihole at the origin: 

*i/3 = n(^-%) 3 exp-Ei^r/K) (6) 

i<j i 

*i/3,,h = n^n(^-^) 3 exp-E|z 4 |V(4/o)- (7) 

k i<j % 

Note that only l enters the wavefunctions and m does not. This is true also of the subsequent 
functions written down by HaldanJ, Halperinii, Jain! and Rezayi- ReaoM 

From these wavefunctions one has gleaned many of the basic features of the FQHE such 
as the incompressibilityS, fractional charge of the quasiparticles§H§'i, and their fractional 
statisticsi'0. These successes notwithstanding, there are many reasons to pursue the hamil- 
tonian approach in parallel. If it can yield these wavefunctions in some approximation, the 
nature of that approximation will be physically revealing. One can calculate other wave- 
functions and correlation functions at general (q, lv), one can couple the system to impurities 
to study the effects of disorder, or to a confining potential to study the physics at the edge. 
The underlying theory will tell us how the charge of the excitations gets renormalized from e 
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in the case of free electrons to e* = e/(2ps + 1) for the principal fractions pj (2ps + 1), and in 
particular how it vanishes for the case of v = 1/2, becoming dipolar 

It should also tell us who carries the Hall current if the quasiparticle (the composite fermion) 
is neutral. It should shed light on the order parameter with algebraic order proposed by 
Girvin and MacDonald^ and studied by ZHK0, and the order parameter with long-range 
order proposed by ReadEl. 

The theory must also address the question of the bare and renormalized masses in some 
detail. The notion that low energy physics can be described with no reference to m or higher 
Landau levels is plagued with subtleties as emphasized by Girvin, MacDonald and Platzman 
(GMP)i. First, the current operator depends explicitly on m while the Hall conductance 
<j xy does not. For m to drop out, one must necessarily invoke virtual transitions to the next 
Landau level@ with matrix elements of order m. Next, one must understand how m gives 
way to m* in the rest of the low energy physics and how the latter is to be calculated in 
terms of the interactions. Even if one begins with a phenomeno logical m*, there are places 
where m enters. For example, the cyclotron pole is determined by m according to Kohn's 
theorem^!. In the Fermi-liquid-like state of u — 1/2, Simon and Halperinil argue that it is 
necessary to bring m back by invoking a suitable Landau parameter f\ and the relation 

m m* Z7T 

a scheme that was briefly alluded to in hlrR Now one may argue that the Kohn mode is 
part of the high energy physics and an effective low energy theory is not obliged to explain 
it. But as Simon, Stern and Halperin (SSH)0 point out, m enters the low energy physics 
as well. Consider the zero point energy of eB/2m per electron, which one may ignore as a 
constant. Suppose we impose on the system a small slowly varying field 5B(r). This amounts 
to imposing a potential eSB(r)/2m, the response to which is surely m dependent. Also, if 
for any reason the density becomes nonuniform, there will be a current due to incomplete 
cancellation of cyclotron currents^: 
g 

jmag = X Vn. (9) 

Both these effects imply that the charge-current correlator K i = K 1Q has a piece 
(iqe/2m)K 00 as q — ► 0. SSH show that these effects can be reproduced by appending to 
each particle a magnetic moment e/2m. It is desirable to have this occur naturally within 
a computational scheme. 

The variant of CS theory we propose here solves many of these problems. 

To set the context, let us begin by recalling the basic idea behind CS theories, using as 
a concrete example v = 1/3 and the goal of getting Laughlin's wavefunction ^1/3. Clearly, 
the more we know of the answer the less we have to calculate. The usual bosonic CS theory 
of Zhang, Hansson and Kivelson (ZHK)i be gins by writing the wavefunction for the state 
as 



#i/ 3 fo 



phase of * 1/3 $>cs(ri) (10) 



and then solving for ^cs which describes CS bosons. At mean field level the bosons condense 
into the zero momentum state, giving ^ C s — 1 an d we obtain just the phase of \J/i/ 3 . Each 
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particle sees three flux tubes in the others, but does not avoid them with the triple zeros, 
a correlation essential to the success of Laughlin's wavefunction. These zeros, which imply 
correlation holes, arise upon going beyond mean-field and including fluctuations^!. This is 
also true for fermionic CS theories of Lopez and FradkirH Kalmeyer and ZhangEI, Halperin, 
Lee and Read (HLR)@ and Kwon, Marston and Houghton0, though wave functions have 
not been explicitly worked out this way for the case of v = 1/2. For the same reason, the 
composite bosons and fermions all carry the full charge of the electrons at mean-field level 
and the fractional charges of the quasiparticles are again to be found in the fluctuations. 
An exception is the theory of Rajaraman and SondhS These authors include all of ^1/3 
as a prefactor in Eqn.([K]) so that at mean field level they obtain all of ^1/3. On the other 
hand, the fluctuations about mean field are described by a complex CS field whose real 
and imaginary parts come from the phase and modulus in the prefactor.Q To go beyond 
mean-field theory using diagrams one must face the fact that the hamiltonian is not her- 
mitian. However the formalism has already been applied successfully to questions involving 
solitonS 

The method we propose manages essentially to append as a prefactor both the modulus 
and phase of ^1/3 (and its generalization to fractions of the form p/{2ps + 1)) without the 
use of complex vector potentials. The correlation holes give the right charge e* to the 
quasiparticles. The theory has two sectors, and with one (desirable) exception, m appears 
in one and m* in the other. 

Setting aside details for later, here is how this is done. Recall (or go and read, if you are 
too young) the work of Bohm and Pinesi! on the electron gas in three dimensions. Besides its 
fermionic quasiparticles and particle-hole excitations, the gas has plasmons in the spectrum. 
These arise as poles in the density-density correlations. The poles have a negligible imaginary 
part (i.e., the plasmons are very long-lived) till the electron- hole continuum runs into the 
plasmon pole at large enough q. Rather than treat the small q plasmons as composites built 
out of the original degrees of freedom, Bohm and Pines elevate them to elementary particles 
with their own set of states and operators in an enlarged Hilbert space. The plasmons 
couple to the electrons in way that reproduces the pole and residue in the density-density 
correlations. Having introduced the plasmons as elementary objects, one must ensure that 
the particles do not duplicate them. To keep the problem unchanged in the enlarged Hilbert 
space, only states obeying a set of constraints (one for each q at which a plasmon was 
introduced) are deemed physical. The constraints freeze all fermionic collective degrees of 
freedom that were assigned to plasmons. Placing the plasma oscillators in their ground 
states gives, upon projection to the physical states, a correlation factor to the fermionic 
wavefunction, which is simply the exponential of the Coulomb potential energy associated 
with each configuration. If the principal coupling between the fermions and plasmons is 
removed by a canonical transformation, it is found that the fermion mass gets renormalized 
upwards while the plasmon pole is fixed at the classical value as q —>■ 0. The number of 
plasmon modes introduced is treated as a variational parameter and is a small fraction of 



1 A CS theory with complex electric and magnetic fields was briefly discussed up by Girvin in the 
closing article in Reference (0). 
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the number of fermions. 

We adapt this scheme to our problem with the magnetoplasmons playing the role of 
plasmons. For a variety of reasons we pick n , the number of oscillators, to equal n, the 
number of particles. Putting the magnetoplasmons in their ground states gives the correla- 
tion zeros of the trial wavef unctions. Our field operators for composite bosons and fermions 
create the electron as well as its correlation hole, the latter being described by collective 
coordinates. When we decouple the oscillators from the particles, 1 / m* gets renormalized to 
zero in the noninteracting limit and picks up a finite value when interactions are turned on. 
The zeroth order approximation to 1/ m* is readily evaluated. The decoupling, performed in 
the long-distance limit and the RPA, also changes the charge of the quasiparticles from e to 
e*. Upon decoupling, the n constraints act on the particles alone, the reduction in degrees 
of freedom being exactly what it takes to get LLL physics. This is confirmed by the fact 
that the density (limited to the sector in which the oscillators are in their ground states) 
now obeys the magnetic translation algebra of GMpi. All dependence on m is contained in 
the oscillator sector with one exception: the particles couple to the external magnetic field 
with a moment e/2m, the need for which was discovered by SSH0. We find that in the 
clean limit, all the Hall current is carried by the oscillators, which takes care of the nagging 
question of who carries the Hall current when the composite fermions become neutral at 
v = 1/2. Our formalism gives a precise meaning to the very useful concept of the composite 
fermion . All this will be explained as we go along. 

The plan of the paper is as follows: In Section II we recall the key ideas behind flux 
attachment in both the wavefunctions and CS approaches. Section III describes our CS 
representation which results from enlarging the Hilbert space and doing a canonical trans- 
formation to reach what we call the middle representation (MR). In Section IV we derive 
correlated ground state wavefunctions for v = l/(2m+l), 1/2 and the Jain series p/(2p + 1). 
(The case p/(2ps + 1) calls for a straighforward generalization and we do not discuss it in 
any detail.) For the first case, described by a bosonic CS theory, we display operators that 
create quasielectrons and quasiholes and a candidate for the Read operators. We show that 
our composite boson and fermion operators create particles in which an electron is bound 
to zeroes, and not just flux quanta. In Section V we perform the canonical transformation 
to the final representation (FR) in which the particles and oscillators are decoupled in a 
long-distance approximation. We exhibit the substantial mass and charge renormalization. 
We pay special attention the gapless case v = 1/2 which arises in the limit p — > oo and 
compare it to the work of HLR. Discussions and conclusions follow in Section VI. 

II. FLUX ATTACHMENT AND CHERN-SIMONS THEORIES 

In two space dimensions one has the remarkable freedom to alter particle statistics by 
attaching point flux tub In particular, a fermion carrying an odd/even number of flux 

quanta will behave like a boson/fermion. This fact is exploited as follows. Let us consider v = 
1/3. In the noninteracting limit, the ground state is macroscopically degenerate. Trading 
the electrons for composite bosons (CBs), each carrying three quanta of statistical flux 
opposed to the external field, one ends up on the average with composite bosons in zero field. 
Alternatively, trading them for composite fermions (CFs) carrying two quanta (opposing the 
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external field) each, one obtains CFs with one quantum each on the average, i.e., enough 
to fill exactly one Landau level. In either case we have a nondegenerate starting point for 
including the effect of interactions. The notion of dealing with the flux on the average to get a 
simple starting point was first exploited by Laughlinll in studying anyon superconductivity. 
Of course one must deal with the (possibly singular) fluctuations to complete the analysis 
of the problem. 

Let us see how flux attachment is implemented in building wavefunctions and CS field 
theories. 

It was noted long ago by Halperini that the 3n zeros of Laughlin's wavefunction are all 
tied to the location of the particles, while only n are required to do so by Fermi statistics. 
The electron and the triple zero are attracted to each other by the Coulomb force, and behave 
like a neutral boson, which condenses in the Laughlin state, as emphasized by ReadH. Thus 
one can write 



*i/3 = Y[(zi - zj) 3 exp 

i<3 



xl-h (ii) 



where xl attaches three flux quanta to the CB whose condensed wavefunction is simply the 
factor of unity. 

On the other hand for fractions of the form pj (2ps + 1) one does not have such explicit 
wavefunctions. Here is where Jain's notion of composite fermion (CF)i comes to the rescue. 
First he notes that we can just as well write 

*i/3 = life - z 3? ex P - E = xl ■ xi- (12) 

i<j i 

The picture now is that the electrons have been traded for CF's which carry two flux quanta 
opposed to the external field, and fill exactly one Landau level in the weakened mean-field. 
The first factor of \\ attaches two flux quanta and the last stands for the filled Landau 
Level of the CF. Now if v = p/(2ps + 1), with p > 1, bosons in zero field are not attainable, 
but Jain's idea composite fermion (CF) idea is still applicable. By attaching 2s units of 
flux opposed to the external one, the electron problem gets mapped on average to that of 
composite fermions which see 

\ = --2s = l/p (13) 

units of flux each, i.e., fill p Landau levels. For example if v = 2/5, we have p = 2, s = 1 
and the Jain wavefunction is 

*2/5 = V X IX2 (14) 

where \2 describes two filled Landau levels and V projects the wavefunction to the LLL. 
(The magnetic lengths appearing in the gaussians in xi an d X2 must be chosen to reflect 
the magnetic fields associated with filling one and two levels respectively.) When v = 1/2, 
(the s — 1, p — > 00 limit) the external and attached fluxes cancel and we end up with 

*i/2 = Vxi \FS) (15) 
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where \FS) stands for the Fermi Sea. This function was proposed independently and studied 
in detail by Rezayi and Read (RR)0. They showed that a gap is not needed for the success 
of the idea. 

It is important to note that there is no unique notion of flux attachment in the wave- 
function approach. For example, instead of using the factor \i t° attach a unit of flux one 
could use a factor that did not contain the gaussian, or even the zero when the particles 
approached each other, by the making the replacement 

(* - *i) - TT^TT (16) 

I i j I 

If the latter choice had been made for u — 1/3 one would have lost two of the three zeros 
that produced the excellent correlations as the particles approached each other. Of course, 
such a wavefunction would have been a much poorer approximation to the true ground state 
than one where the correlation zeros were present. There is no reason to make the poor 
choice in the quest for wave functions, but as we shall now see, the choices are rather limited 
in the CS approach to which we now turn. 

In contrast to the, flux attachment in CS field theories is a lot more restrictive, and is 
usually carried out as follows. In first quantization one introduces a wavefunction for the 
CS particles in terms of the electronic wavefunction as follows!! 

*e = II JTZTTJ^os. (17) 

i<j K« 4, 3\ 

where I is the number of flux quanta to be attached. The CS particles are to be quantized as 
fermions/bosons, for I =even/odd since the prefactor produces a factor (—1)' under particle 
exchange. The prefactor introduces a gauge field a in the Schrodinger equation for ^cs- 
The field obeys 

is 5 -?<■('-«>■ < 18 ) 

In second quantization the CS transformation is represented by the operator relation 

z x (r — r 



V> e (r) = exp 



fc(r') (19) 



where ip e and ipcs are the field operators that annihilate the electron and the CS composite 
particle respectively. (They are not to be confused with the wavefunctions $). Also, 

p(r) = V!(r)V e (r) = ^ cs (r)^ cs (r) (20) 

is the density. The hamiltonian(density), before turning on electron-electron interactions 
(which will be added later), is 

t |HV + e A + a)| 2 

Hcs = r C s 2^ ^ cs ( > 
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where A is the external vector potential, m is the bare mass and 



V x a + . , 

— = ,H: (22) 

Since A describes a magnetic field pointing down the z-axis, and the CS field has its 
flux of I quanta per particle pointing up, they can be played off against each other. There 
are two kinds of v for which things are particularly simple. If v = l/(2s + 1), by attaching 
2s + 1 quanta per particle, we can trade the electrons for bosons that see zero net field 
on the average. If v = l/2s, by adding 2s quanta per particle we can trade the electrons 
for fermions in zero average field. Hereafter we use the convention that the average values 
have been removed from a and = p (the density) both of which henceforth refer to 
normal-ordered quantities. (The normal ordering will occasionally be made explicit). The 
bosons condense in zero field and the fermions form a Fermi sea (FS), both of which are 
nice nondegenerate starting points explored by Zhang et a (bosonic) and Kalmeyer and 
ZhangS, HLR0, Kwon et a£3 (fermionic). At the mean- field level the wavefunctions one 
obtains are 

*vi = n TT^TF*^ (23) 

i<j 

where ^cs is the Fermi Sea (FS) for v — 1/2 and ^cs — 1 for v = l/(2s + 1). The multiple 
zeros expected in the correlated wavefunctions are absent. To obtain these one needs to 
analyze fluctuations about the mean fielcM 

When v — pj {2ps + 1), by attaching 2s opposing flux quanta per electron, we can cancel 
enough of the external flux so that what remains, eA*, is enough to fill exactly p Landau 
levels : 

eB 

eB* = — — . (24) 
2ps + 1 y ' 

These are the fractions considered by Jain!. Indeed Jain's case reduces to the two 
previous ones when p = 1 and o = oo respectively. The corresponding field theory was 
analyzed by Fradkin and LopezM who computed some response functions and also showed 
that fluctuations produced the modulus of the Laughlin function in the limit of particle 
separation large compared to l . 

Let us summarize: In the quest for wavefunctions we can build in not just the phase, but 
all of Yl(zi — Zj) 1 , or even the ubiquitous gaussian factors into the process of flux attachment. 
In contrast, the standard way of attaching flux in CS theories leads only to the phase of the 
zeros at mean-field level, while the zeros themselves emerge from the fiuctuationsEl'0. It is 
possible to obtain the zeros in a field theory at mean-field level: however, doing so leads to 
a complex vector potential and a non-hermitian hamiltonian, as pointed out by Rajaraman 
and Sondhi who analyzed this possibility in some depth!!. Let us now turn to our version 
of the CS theory and see how correlated wave functions are obtained. 
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III. OUR CS THEORY 



Let us recall the CS hamiltonian (density) before introducing any interactions: 

^ = ^ IHV+ 2 ^ + a)|2 fe. (25) 

Here a refers to the fluctuations about the mean (which has been used to reduce A down 
to A*), and obeys 

V x a . „ . 

— =:*V:. (26) 

To make explicit the fact that a is really a dependent field, let us rewrite this as 

H = -^fcs(-W + eA* + (Vxj-^ri^fc (27) 

where the inverse curl is uniquely defined by its Fourier transform if we demand that the 
answer be transverse. If you find this form of H forbidding, it is our intention; Eqn. (p25|) 
hides the complexity of the CS hamiltonian. 

We are now ready to enlarge the Hilbert space. Consider a disk (in momentum space) 
of radius Q = hp, the Fermi momentum of the electrons, which contains as many points as 
there are electrons. For each q in this disc, we associate a canonical pair of fields a(q), -P(q) 

[a(q),P(q')] = (27r) 2 5 2 (q + q'). (28) 

These define a pair of longitudinal and transverse vector fields 

P(q) = zq P(q) a(q) = — iz x q a(q). (29) 

The Hamiltonian density of Eqn. fl2"T| ) is completely equivalent to 

H = 7T-ftcs(-W + eA* + a + (Vx)-^^) 2 ^ (30) 
provided we restrict our selves to states in the larger space obeying 

a(q)\physical) = q < Q. (31) 

In other words, [H, a] = allows us to find simultaneous eigenstates of H and a, and the 
eigenstates with a = solve the original problem. 

Let us note for future use how we are to extract wavefunctions for the CS particles from 
any solution in the bigger Hilbert space. Let ^csi x ) denote the CS wavefunction. In the 
larger space we have, in obvious schematic notation, the following resolution of the identity 

dxda\x a) (a x\ (32) 

where x denotes all the particle coordinates and a stands for a(q), [0 < q < Q]. The 
projection operator to the physical sector is 

p — J dxda\xa) (ax\-^-^. (33) 
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This means that if *&(x, a) is a generic wavefunction, then the projected version is 

a) = a) = 0). (34) 

Indeed one can show that every physical vector must have such a form. 
Let *&p(x, a) obey the eigenvalue equation (in first quantization) 

H(x,p,a) ^p(x, a) —E ^/ p (x,a) (35) 

where H(x,p,a) is the hamiltonian with a added on as in Eqn.flHUp. Such eigenfunctions 
will exist since H commutes with the constraint. Feeding in Eqn.fl34]) for the wavefunction, 
moving the delta function through H to its left using the commutativity, and canceling it 
from both sides, we obtain 

H(x,p, 0) ty(x, 0) = E ty(x, 0) (36) 

establishing that 

^cs(x) = V(x,0). (37) 
In other words "&cs is embedded within X) ft) clS follows: 

*p(x,a) = 5 M*cs(x). (38) 

Let us now return to the Hamiltonian density of Eqn. ( j30|) which contains the very nasty 
inverse curl. We deal with it through the following unitary transformation 

U = ex V [J2iP(-q)—p(q)}. (39) 

q Q 

In the above and what follows q sums stand for integrals 

and the vector nature of q is often suppressed. Under the action of U 

rfa{q)U = a(q) - (41) 

Q 

Q 2tt1 

ipcs(x) = TPcp{x)e X p[J2iP(-q)—e- iqx } (42) 

q Q 

^ cs {x){-iV)^ CS {x) = ^ CP {x){-iV + 2vr/P(x)) ^ CP {x) (43) 
H = —^ CP {-iV + eA* + a + 2nlP + dafipcp (44) 

= (a - ^^-)\physical ) 0<q<Q (45) 
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where Sa refers to the dependent short range vector potential (Vx) _1 27r/p for q > Q that did 
not get cancelled by the unitary transformation, and ipcp refers to the composite particle 
(boson or fermion) field. The above description of the problem, after the canonical trans- 
formation by U, will be referred to as the middle representation (MR). We argue that ^cp 
better describes the composite bosons and fermions alluded to in the literature than does 
^cs- Whereas the latter are associated with particles carrying just flux tubes, the former 
are associated with particles that carry flux tubes and the correlation holes, i.e., describe 
electrons bound to zeros, as we shall see. 

Suppose we have a physical wavefunction ty*f R (x, a) in the MR: 

*? R (x,a) = 6{a ~^ ) * MR (x,a) (46) 

What is the corresponding CS wavefunction? Given that U implements a translation of 
a(q) by it follows that 

»_(*, a, = *r («. - + ^> = H* M *(*, ° + 2j v ] = W)* MR(X ' 2j t ] (47) 

from which follows, upon invoking Eqn.(^), the invaluable result 

*cs(*) = * M V— )■ (48) 
Q 

The hamiltonian in Eqn.(fl4T) 

H = —ib P {-iV + eA* + a + 2tt/P + Sa) 2 ip C p (49) 
has a local gauge invariance under time-independent transformations: 

^cp - e 2 ^ CP (50) 
P -> P-VA (51) 

where A has only Fourier modes with q < Q. The constraint merely states that physical 
states must be singlets of the generator of these transformations: 

&^-:^ip:)\physical)=0. (52) 

We may also write the above as 
qa 



2tt/ 



: p(q) :)\physical ) = < q < Q. (53) 



Had we chosen Q to be infinite, <5a would have vanished and we would have entered what 
is called the Weyl or a = gauge. In the Bohm-Pines case, Q is very small compared to 
kp and chosen as a variational parameter^!. In our problem the choice Q = hp recommends 
itself over all others repeatedly, as we shall see. The gauge we use is a nonstandard one, tailor 
made for this problem. Sometimes one is asked how a choice of gauge could matter, given 
that the physics is gauge invariant. This misses the point, which is that one can use the gauge 
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freedom to highlight some specific feature of the theory. Thus in Yang-Mills theories, the 
Coulomb and unitary gauges display the physical degrees of freedom, the Feynman-Lorentz 
gauges the Lorentz invariance, the gauge the renormalizability etc. 

We now explore our hamiltonian Eqn.(f44"D, expanding it as follows (and dropping the 
subscript CP): 

H = — |(-iV + eA* + a + 27r/P + 5a)^| 2 (54) 

= 7 ^|HV + eA>| 2 + ^(a 2 + 4vr 2 / 2 P 2 ) 
2m 2m 

+ (a + 2vr/P) ^U-iV+eA*U 

2m 

+ i|^(a + 2 7 r/P) 2 
2m 

+ 7~V> + H V+eA*)^ + ^(2(a+27rZP)+5a)-5a (55) 

2,771 2j 777 

= Ho + Hj + Hn + Hsr (56) 

in obvious notation, H sr being the terms associated with the nondynamical short-range 
gauge field 5sl. Note that we are yet to add interactions. 



IV. WAVEFUNCTIONS FOR GROUND STATES AND EXCITATIONS 

In this section we will show how the above hamiltonian leads to some well known cor- 
related wavefunctions upon making the simplest approximation. Let us first consider the 
fractions v = l/(2s + 1) and l/2s for which A* = 0. After this we will turn to Jain's 
principal fractions p/(2ps + 1). We will focus on the case s = 1 since larger values of s do 
not tell us anything new. 



A. v = 1/3 

For this case we trade the fermions for bosons by attaching I = 3 flux quanta. Let us 
focus our attention on just H in Eqn.ffib]). It reads 

^ = 2^ IHV) ^ |2 + 2^ (a2 + (67rP)2) - (57) 

There are no cross terms between a and P since they are transverse and longitudinal 
respectively. Given that (a , P) are canonically conjugate, we see that the second term 
describes harmonic oscillators of cyclotron frequency 

67m eB 

= — = 0Jc- (58) 

m m 

The oscillator hamiltonian may be written as 
Q 

H osc = Y d A\q)A{q)u c (59) 
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where 



A(q) = (a(q) + QiriP(q)) j^fVhx. (60) 

We are now ready to discuss the upper limit Q that controls the number of plasma 
oscillators. Given that there are only n electrons in a plane with a total of 2n degrees of 
freedom, the number of independent oscillators is bounded by 2n. We choose the number of 
oscillators to be introduced to be equal to the number of electrons, i.e., Q = k F . There are 
many reasons for this choice, only one of which is intelligible at this point. Recall that we 
trying to describe LLL physics for the low energy sector along with high-energy cyclotronic 
physics. The oscillators clearly correspond to the latter. To pay for them, the particle sector 
has to give up n = n degrees of freedom, which is precisely what it takes to project from 
the full Hilbert space to the LLL. Although mathematically any value of Q is allowed (being 
different choices of generalized gauges) our choice will prove the most suitable. 

Our quest for the ground state of the full interacting problem can be described in two 
equivalent ways. 

Firstly, we can try a variational approach using a simple product wavefunction in the 
enlarged Hilbert space, initially paying no attention to the constraint. The reason is the 
following. Suppose one were capable of finding the average of the Hamiltonian in arbitrar- 
ily complicated variational wave functions. Then the one with the lowest energy in the 
enlarged Hilbert space will be the true ground state. However, since the Hamiltonian is 
gauge invariant, we expect the ground state to be gauge invariant as well, since we know 
of no case where gauge invariance is spontaneously broken. Thus in the variational quest 
for the ground state, we can forget gauge invaraince; the winner will have that feature. Of 
course, we can can only carry out the variational procedure for product states, so all the 
above considerations can be expected to hold approximately. We will use for this purpose 
the ground state of H which we can write down by inspection: the composite bosons all 
condense into the zero momentum state with wavefunction ^/ C b( x ) = 1 and the oscillators 
occupy their ground states: 



^ osc = exp 
giving us the MR wavefunction 



(61) 



® MR (x, a) = exp 



(62) 



Not only does this minimize Hq, both Hi and Hn have zero average in this state: Thus 
this wave function minimizes H + Hj + Hjj among product wavef unctions. In addition, any 
density-density interaction (which will go as a 2 q 2 V(q) since p(q) ~ qa by the constraint) 
and H sr will make some difference only at very short distances, where the wavefunction will 
be found to have excellent correlations that keep the particles apart. 

Equivalently, we can think of the product wavefunction as a zeroth order starting point 
for a perturbative solution based on H . From this vantage, the approximation has a chance 
of being relevant only in the presence of electron-electron interactions: without them there 
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is no unique ground state in the exact solution and the product ground state of H must 
necessarily get destabilized by the neglected terms. Thus interactions are needed for our 
approach to work, even though they are not being explicitly treated. Note also that H does 
not correspond to free electrons; for that we need all of H in Eqn. (|56f) . By dropping all 
terms but Hq, we have gone from the free problem with degeneracies to one with a unique 
ground state. 

There is just one problem: the wave function in Eqn. (|5^) is in the enlarged Hilbert space, 
whereas the physical wavefunction must be constrained to lie in the physical subspace. We 
implement the constraint by simply multiplying the above wavefunction by a delta function 
of the constraint, i.e., by projection. As explained in the steps leading up to Eqn.(f48|), ty C s 
is obtained by setting a(q) = 6^p(q)/q in ty MR (x, a): 



cs 



exp 



exp 



exp 



^3tt : p(q) : — : p(-q) : 



^ i 3 

2n'^2 J dx\n\x — Xi\+ constants) 



exp 



i<j 



(63) 
(64) 
(65) 
(66) 



The steps connecting to the penultimate and ultimate lines are from Kane et «fl. The 
integral over x in the penultimate line may be interpreted as the potential energy of a 
point charge at Xi, due to a uniform charge density Aim and the two-dimensional Coulomb 
potential V(\x — y\) = In \x — y\. (To find the potential at x, it helps to first find the 
field by invoking the two-dimensional Gauss's law § E • dr = J d 2 xp(x) and then integrate.) 
Putting back the phase factors from Eqn.([T7|), gives us Laughlin's ^1/3. 

Let us recapitulate. We began by showing that given an exact gauge invariant eigen- 
function of the enlarged hamiltonian, we could get the physical one by dropping the delta 
function of the constraint that it must necessarily have as a prefactor ( Eqns. (|38D and 
(48|). ) In the end we found an approximate variational wavefunction fy MR (x, a 



w ! '' ' - - A that one 

could argue had a good energy, but was not gauge invariant, i.e., did not vanish outside the 
constraint. We then projected it to the physical subspace by appending the delta function 
${x) (x being the constraint) and evaluating the product wavefunction on the constrained 
subspace to get ^cs- Could not multiplication by S(x) ruin the good energetics? In general 
yes, but not here: given that ^f MR (x, a) is an approximate eigenstate of the exact H, i.e., 
obeys the wave equation 



MR 1 



EV MR (x, a) 



(67) 



2 Due to the cut-off Q, the form of \& we get is good only for \z% — Zj\ » 1/Q- To continue this 
form all the way in to get a triple zero, we must invoke the LLL condition. 
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we can premultiply both sides by <5(x(g)), commute it through H to verify that 
6(x(q))^ MR ( x i a ) obeys the same equation. 

Our trick for extracting the ground state can be extended to find a few other states 
as follows. Consider the sector in which there is a vortex at the origin. Among all such 
electronic wavefunctions there must be one with lowest energy, the ground state in this 
sector. Let us then write electronic wavefunction as 



q, vortex _ TT e »E/i^ 



(68) 



and subject ^f' e , which is free of this vortex, to the previous CS transformation. Instead of 
Eqn. (|30|) we will end up with 



H = J_^t(_iV + eA* + a + (Vx^Gnp + — 
2m r 



(69) 



the last term in brackets being due to the vortex prefactor. 

The unitary transformation, which must now get rid of the last two terms in the bracket, 



is 



U = exp[J2iP(-q)(-p(q) + -)}. 

q Q Q 



(70) 



The bosonic wavefunction will again be unity since the external flux is zero on average, 
including the point flux tube at the origin. The oscillator wavefunction will still be the 
gaussian, 



exp 



1 



-E — « 2 (q) 

^ 12tt VH; 



but the constraint will now be modified to 



a(q) 



67rp(g) 2tc 
q q ' 



(71) 



(72) 



because of the modified U. 

It is clear that we have now placed a phantom charge of size 1/3 at the origin. The 
plasma will screen it and produce a correlation hole of charge —1/3. This is verified by 
evaluating the corresponding wavefunction using the same procedure as before to obtain 



r)\e 



Zj I exp 



i<j 



EN7^ 



(73) 



where we have displayed the answer for the case in which the vortex is at rj rather than the 
origin, to show how even the gaussian factor associated with its location appears naturally. 
Note that the magnetic length associated with it obeys Iq 2 = 3l 2 , appropriate that of a 
charge —1/3 object. Going back to the electronic wavefunction, \zj — r]\ will become Zj — t] 
giving us Laughlin's quasihole state, again with the understanding that we continue down 
to z — rj = using the LLL condition. 
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To create two quasiholes, at 771 and 772, we do the obvious extension of the preceding and 
find 



11/3 p -(\m\ 2 +\v2\ 2 )/uf 



n \ z j-vi\ n \ z j~v2\ ni^ _z ii 3ex p 



i<j 



Observe that the wavefunction contains the normalization factor 1 771 — 772 1 1//3 of the two- 
quasihole state. If we follow Halperini and drop the mod sign (by a singular gauge trans- 
formation) we can get a pseudo wavefunction for the quasiholes and infer their fractional 
statistics. 

It is possible to assign operators that create the quasiholes in the following sense. Con- 
sider the action of the operator 



ip(0) = exp 



q y 



on the unprojected ground state which we denote by 

\Osc = 0,^ C B = l). 
Its action is one of translation in a: 

e -a 2 (g)/127r _^ e -ia(q)+2f] 2 /127r 

If we now obtain the projection, pip(r])\Osc = 0, ^cb = 1): we end up with 



cs\ x > Q ) = IIN - Zj\ 3 exp 

3 i<j 



£N7^ 



(75) 



(76) 



(77) 



(78) 



It follows that the operator which creates a quasihole at the point 77 is 



) Z S 1 



exp 



2wi 



E—P(-1> 



,iqrj 



i q 



Q 



(79) 



where 8j tV is the angle between the x-axis and the vector connecting rj to Zj. 

Note that the projection is done only at the very end. In other words, to create two 
quasiholes, we need 



Pipqhfa) fpqh(Vl) \0SC = 0, ^ C B = 1) 

and not 

p4>qh{m)p 4>qh(vi) \Osc = 0, ^CB = i) 



(80) 



(81) 



In fact p^pqhp has zero matrix elements because ip q h is not gauge invariant. 

The adjoint operator V^l 7 ?) will append to the ground state wavefunction the factor 
Ylj(zj — v) 1 which will indeed have a charge 1/3. However the resulting wavefunction does 
not belong to the LLL since negative powers of z are not allowed. Thus our formalism does 
not solve this problem with quasiparticles, first encountered by Laughlin. 
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Let us summarize what has just been done. In the CS representation, the hamiltonian 
did not contain any spurious degrees of freedom, but was very nonlinear due to the (Vx) _1 
terms that corresponded to the nondynamical field a. We moved to the MR in which 
the oscillator degrees of freedom were introduced. The hamiltonian, which now contained 
new canonical degrees of freedom, but no long-distance (q < Q) nondynamical fields, had 
invariance under time-independent gauge transformations. Physical states had to be gauge 
singlets to compensate for the new degrees of freedom. We found an approximate variational 
product ground state \Osc = 0, ^cb = 1), which was the ground state of H , a nongauge 
invariant part of H with a unique ground state that could be read off by inspection. We then 
projected this ground state using p to get Laughlin's wavefunction tyy 3 (x). This procedure 
succeeds only because of particle-particle interactions which render the true ground state 
unique and nondegenerate. Continuing, we saw that ip q h acting on \Osc = Q,^cb = 1) 
produced states, which, upon projection, led to normalized quasiparticle states. It appears 
that there is a dictionary between the eigenstates of H and the full interacting hamiltonian: 

\Osc = 0,^ C B = l) ^ *i/ 3 (ar) (82) 

4>qh(llk) ■ ■ ■Tpqh(?ll)\C>SC= 0,^ C B = 1) ^1/3, qh(x,Vl,- ■ 'Ilk) (83) 

All these states may be viewed as ground states in sectors with a given number of vortices. 
We do not know if the correspondence carries over to excited states. 

We can now understand how we have managed to get the phase and modulus of ^1/3 in 
the middle representation. Recall Eqn. (f42|) 

4*0*0 = 4s(z)exp[£zP(-?)-e-n (84) 

q 1 

Thus ipCB creates a CS boson and introduces a triple quasihole of charge —1, thereby 
creating the extra particle along with its correlation hole. Starting from Halperin's original 
observation, all physical pictures, particularly Read'sBS, have emphasized that electrons 
bind to zeros of the wavefunction, not just the vortices in the phase, i.e, the flux tubes. In fact 
only a zero with an accompanying charge deficit will attract an electron. This is the reason 
we prefer to use the term Composite Boson for the particle created by ipcs rather than by 
ipcs- The fi rs ^ f ac t° r °f V'aB carr i es the flux, i.e., the phase, and the second, the magnitude, 
i.e., the correlation zeros, of ^1/3. (The latter has the effect of introducing the hole only 
upon projection.) Crucial to the success of all of this is the assignment of independent 
degrees of freedom to collective coordinates, rather than trying to express them in terms of 
particle coordinates. It is only after a(q) became a canonical coordinate whose conjugate 
momentum P could be used to translate it, that we could so readily produce the correlation 
hole. 

The preceding arguments suggest that is a candidate for being Read's operator!!]. 
Let recall his reasoning for the choice of condensate. If we take a Laughlin state with N +1 
electrons and kill one at z, it looks like an N particle state with a unit charge correlation 
hole there. Schematically 

(N + l\4 lectron (z)U(z) 3 \N)^0 (85) 
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where U(z) in first quantization appends a factor Yii( z ~ z i), i- e -> creates a quasihole at 
z. Thus 4>ii ec t ron (z)U (z) 3 condenses. Given the phase relationship between ipeieictron and 
ipcs-i this means il)\j S {z)\U{z)\ 3 condenses. Comparing to Eqn. (|84|) suggests ipcB is the Read 
operator. There is, however, the caveat that it acts in the enlarged Hilbert space rather 
than the physical subspace. While it is true that 

pexpEiP(-g)— e- iqz ]\Osc = 0, = 1} = J[ \z - % | 3 * 1/3 = \U(z))\ 3 V l/3 (86) 

q Q j 

it is not true that exp\Y^ iP(— q)^-e~ iqz ] acting directly on the physical state ^1/3 produces 
the correlation hole. Thus it is the preimage of \U\ 3 (in the big space, before projection), 
but not an operator which acting on ^1/3 produces this factor. 

There is another related way to say this. As far as Hq, which is a nongauge hamiltonian 
goes, it is indeed true that 

(Osc = 0,^ C B = 1\^ cb {v2)^cb{Vi)\Osc = 0,^ C B = 1) -> 1 as \ Vl - r] 2 \ ^ 00. (87) 

On the other hand, within the gauge theory defined by H, ipcB is a gauge dependent operator 
which cannot have a mean value in the true gauge invariant ground state. (Introducing the 
line integral between the points 771 and 772 does not help; it m akes it the correlator of ipcs-, 
which has the algebraic order found by Girvin and MacDonaldiS. ) 



For v — 2 we choose I = 2 flux quanta, cancel the external field, and end up with 
fermions in zero field on the average. An analysis of Hq as in the previous subsection will 
give us (upon putting back the phase factors of the CS transformation) 



K3 



\FS) 



(88) 



which differs from the Rezayi- ReadS wavefunction by the absence of projection to the LLL. 
This issue will be addressed in the next section. 



C. The Jain series v 



2p+l 

In this case we we follow Jain! and attach two units of flux per particle to reduce eA 
down to eA*. Since A corresponded to \ jv flux quanta per particle and A* to l/p quanta 
per particle 



A* 1 



A 2p+l 
Now Ho takes the form 



(89) 



S T) 

H = ^IHV + eA>| 2 + —(a 2 + (4ttP) 2 ). (90) 
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First consider the oscillators. They are not at the right (cyclotron) frequency but at 



Aim 



cu = 



2p 



m 



2p + l 



-UJ r = 2uUJ r . 



(91) 



Nonetheless they lead to good wavef unctions. Their ground state wavefunction, upon pro- 
jection, gives a factor 



Yl \zi - Zj\ 2 exp 

i<j 



-E 



2u\zj\ 



(92) 



The derivation proceeds exactly as in the case of v — 1/2, except that we replace Ann by 
2v/ll 

The particles fill p Landau levels. Let us write their wavefunction as 



*cf = X P = exp 



-E 



All 



*2 



f P {z,z*) 



(93) 



where / is a polynomial and l*, 2 = (2p + 1)1% . This leads to a wavefunction 



^cs = II \ z i ~ %| 2ex P 

i<j 



E 



r 4/ 2 



f P (z,z*). 



(94) 



It is very gratifying that the two gaussian factors combine very naturally to give the right 
magnetic length because of the relation 



2z/ + 



2p+l 



1. 



(95) 



Note that there is no further projection to the LLL. (Jain has pointed out that even 
without the projection, the wavefunction is primarily in the LLL in the sense that the 
kinetic energy per particle is only slighly higher than %uo c . ) In any event we consider the 
question of projection to the LLL next. 



V. TRANSFORMATION TO THE FINAL REPRESENTATION 

The MR has yielded many wavefunctions with good correlations. There are, however, 
some weaknesses. First of all, the particles have too much kinetic energy: in the absence of 
Coulomb interactions, they should have no kinetic energy: 1/m* = 0. There is as yet no 
evidence of this mass renormalization. The oscillators, on the other hand, have the right 
frequency only when A* = 0, i.e., for v = l/(2s + 1), 1/2. (Despite all this, the wavefunctions 
are good because the particle mass and oscillator frequency do not enter them, only the 
magnetic length does. As a result, as long as there is some nonzero 1/m* in the end, these 
wavefunctions remain good. In the absence of interactions, 1/m* indeed vanishes and the 
nondegenerate ground state of Hq becomes irrelevant. However, with interactions, it can be 
a good approximation to the true ground state). 
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Another problem has to do with the charge of the particles. The composite fermions are 
supposed to have a charge e* = e/(2p + 1) for v = p/{2p + One way to obtain this 

result is to imagine injecting a composite fermion by adding an electron and adiabatically 
injecting the two quanta of flux. During the injection, there is tangential electric field. While 
the physics may be complicated near the fermion, on a circle of large radius the effect of this 
field is to produce a radial current determined by the bulk Hall conductance. The integrated 
charge outflow to infinity or to the edge, i.e., the charge of the correlation hole, is readily 
found to be — 2ve. Adding to the charge of the injected electron, the composite fermion 
charge becomes 

e * = e (l - _^_) = — — . (96) 
v 2p+l ; 2p+l K J 

The same argument predicts that the composite boson in v — 1/3 is neutral since in that 
case we add three flux quanta and 3z/ electronic units are driven to infinity. 

For v — 1/2 these arguments imply neutral composite fermions. Readi-i gives more 
specific details based on the analysis of 



i^RR = V J[(zi - Zj) 2 exp 

i<j 



\FS). (97) 



where \FS) denotes the Fermi sea. Without the \FS) factor, each electron sits in the middle 
of the double correlation hole and the complex is neutral. The role of the Slater determinant 
is to produce the antisymmetry. It contains factors of the form e* k r = e l ( kz *+ k * z )/ 2 _ Since 



z* acts like 2ZqJ^ upon projection by 7^, its action on the Jastrow factor is to split the 
double zero and move the electron off the center of the correlation hole. The shift in z is 
ikl\ and the dipole moment associated with a fermion of momentum k is el\ z x k. Note 
that projection to the LLL is essential to the above picture. 

The reason all these features are not evident in the MR is that the particles are coupled 
to the oscillators and neither is a true quasiparticle with its ultimate (or nearly ultimate) 
mass, charge etc. One way to deal with this is to integrate out one in favor of the other. We 
take the more symmetric route of eliminating the coupling between them.This will be done 
within the following approximation scheme: 

• We will work at long distances. Thus if any quantity has an expansion in powers of q 
we will keep just the leading term. This involves the neglect of derivative couplings. 

• When we encounter the density operator in a product with other operators, we will 
use the RPA: 

e^- fc »^ ~ n(2n) 2 5 2 (q - k). (98) 



Thus we will ignore H sr , the short range piece due to 5a, and the non-RPA term Hjj. 

The £ 
fermions. 



The analysis of v = ^-j- will be carried out in first quantization. The particles are 
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The hamiltonian H + Hj reads 



# = E^ + E A] (Q)Mq)"o + 9w [J(q)A(q) + A^q)c(q) 



j 2m i 



= T + H osc + Hi 



where 



n 

V x eA* 
A(q) 

e 

c(q) 
V + 



p + eA 

-eB* = 
1 



eB 



8tt 



2p + 1 
[a(q)+AmP(q)] 



Ann 
Aim 

m 



2p 



2p+l 



-iqxj 



A(q),A\q' 
It is useful to know that 
[U_,U + } = -2eB*, 
and that in our approximation 



V X i iVy 

(27r)V(q-q') 



c(g),c T (g) 



-i(q-q')xj 



+ 0(q) 



. gU + E[n J + ,n! 

~ 2e5*n(27r) 2 5 2 (q-q')- 
The other commutators vanish to this order in q: 

[c(q),c(q')}=[J(q),J(q')]=0(q,q'). 
The coupling Hj will be eliminated by a canonical transformation 



U(X ) = e iSoXo = exp 



Q 



\ 9Y,(c\q)A(q) - A\q)c(q)) 



(99) 

(100) 
(101) 

(102) 

(103) 

(104) 

(105) 

(106) 
(107) 

(108) 

(109) 
(110) 

(111) 
(112) 



where A is to be chosen appropriately. The operators of the MR, called Vl old , will be related 
to those of the FR (with no superscripts) by 



£} old = e~ iS ° Xo £le i ^ oXo 
We will also define 
Q(A) = e - iSoX Qe iSoX 



(113) 



(114) 
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so that 

n(Ao) = fi oM n(o) = n 

We start with the flow equations 
<L4(g,A) 



dA 
G?c(g, A) 

dA 



-0c(g,A) 
2eB*n6A(q,\). 



(115) 

(116) 
(117) 



It is understood that < q < Q in the above and what follows. The flow equations have the 
solution 



A(q, A) = cos/iA A(g) sin/iA c(q) 

H 

c(g, A) = cos/iA c(g) + ^ sin/iA A(g) 



/i- 



2eB*n6 2 



1 

2p' 



Consider the kinetic energy, which may be written, using Eqn. 

v ninj eg* 

^ 2m ^ 



as 



2m 



(118) 
(119) 
(120) 

(121) 



The second term, which does not evolve to the order we are working in, will turn out to 
describe the magnetic moment coupling of SSH@. The first evolves as per 



dT eB*9 



dX 
so that 



m 



^(At(g,A)c(g,A)+ct(g,A)A(g,A)) 



(122) 



T(A)=T + ^f: 



rn 



sin2A/i 
2/i ' 



i (q)A(q) + A^q)c(q)) + 



1 — cos 2/iA ,/i 



e 



(^(q)A(q) - -ct(g)c(g)) 

2/i V jJL 



(123) 



The oscillator hamiltonian and Hi assume the following form in terms of FR operators: 



Q / n q2 \ 

H osc = u o yi AUq)A(q) cos 2 a\ - (AUq)c(q) + cUq)A(q))— sin 2/iA + — sin 2 fj,\cUq)c(q) (124) 



and 



Hi U 

cos2/xA(c t (g)A(g) + A\q)c{q)) + sin 2/iA (^A f (g)A(g) - -c t (g)c(g)) 



■ (125) 



We will now choose A = A such that cross terms between the particles and oscillators vanish. 
Gathering the coefficients of A\q) c(q) from T, H osc and Hj we find that they add to zero if 



tan A /i = /i 



1 



2p 



(126) 
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This completely fixes the canonical transformation. 

We now ask what the frequency of the oscillators is and find that the coefficient of 
the A^A term is exactly oo c = eB/m\ This agreement with Kohn's theorem confirms the 
soundness of our decoupling transformation whose only free parameter has already been 
chosen. 

What about the particles? The c*c term has coefficient 
kinetic energy is 

TT- 7 T\ j 



2m 



E 



-VVVn i e" i,(lrIi) n' 

™„ - 



2m 2mn 

3 %3l 



Thus the total particle 
(127) 



We also need to transform the constraint 
qa old (q) 



p° ld {q) 



47T 



on physical states by transforming both sides of the equation. The flow equation 
d ^ = -^=(A( q ,X )+ AH- q ,X)) 



(128) 



(129) 



can be integrated to give 

V87T \ n 

while previous results for A tell us that 



(A(q) + A\-q)) - 



cos fJ,X )(c(q) + c\-q)) (130) 



qa old (q) 
An 







cosfiX (A(q) + A^(-q)) sin/iA (c(g) + c'(-q)) 



;i3i) 



Remarkably, the terms involving the oscillators match on both sides upon using tan fj,X = \i 
and the constraint involves only the particles: 



p(q) 



Uq(1 — cos //Ac 



cos /jlXq sin /xAq ~" 



J> x n,) e -^ = 



ill 



cos /iAo(l + cos /iAo) ~" 



J> x n^e"^ (132) 



This is very fortunate since it is no use decoupling the particles in the hamiltonian if the 
constraint still couples them. 

Now we calculate the transformation of the current 



n 



Jf(q) = V tli e -ii*i + -^q + A(q) 



-, old 



(133) 



to the FR. (Note that the non-RPA part of it has been dropped). We find 



2% cos fj,X 



A(q) = C ^ X ° A(g). 



2vr 



(134) 



Remarkably the entire contribution comes from the oscillators; the part proportional to 
II just cancels out. 
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To summarize, in the FR, for the fraction v = p/(2p + 1) 



ff-E^+E^-^EEEnie-^->n i+ i:^( 5 )^ 5 K(i35) 

il 2 

X(q) = p{q) + ttt^ rr £(? x n^e"^ = (constrain*) (136) 

COS //Aq(1 + COS flAo) j 

jf(q) = J. +W ° A(q) = (J M (-q)y (137) 
y/Iir cos yUA 

p0W(g) = 9£2^o + A t(-g)) + p(g) - - ^ g(g x ILje"^) (138) 

V o7T 1 + COS flAQ 



J 

tan/xAo = = /i (139) 

v 2p 



where once again we are assuming < q < Q. Outside this region there is no difference 
between old and new variables. Note that the transformations do not require that Q have 
any particular value0. 

Consider Eqn.( |i~35D for the hamiltonian, focusing on the particle sector. We see that 
decoupling the oscillator has lead to the third term. In it is a one-particle piece corresponding 
to i = j, those with i ^ j being additional interactions. When we combine this i = j term 
with the first term, we see the mass gets renormalized to : 

1 1 1 Q 

- = -(l--V). (140) 

m* m n q 

As the upper limit Q grows, the particle mass increases and the kinetic energy gets 
quenched. In the noninteracting case we are studying, the kinetic energy must be fully 
extinguished. Note that our choice Q = hp, i-e., n osc = n accomplishes this exactly. In 
addition, this is the right number of constraints, all in the particle sector, to restrict it to 
the LLL. A smaller choice of Q gives the particles kinetic energy of order 1/m and also does 
not impose enough constraints to project to the LLL. One is free to use a smaller value for 
Q, but one will have a harder time going from such a hamiltonian to the correct behavior 
for the noninteracting case. A larger value Q > kp leads to a negative 1/m*, clearly a 
bad starting point when we go on to introduce interactions. Note that the kinetic energy 
of the fermions is quenched independent of their momentum. We do not know what set 
of diagrams, if any, corresponds to carrying out this canonical transformation. Recently 
Chari, Haldane and Yangpil worked on a gauge invariant approximation scheme for mass 
renormalization in a theory with flux tubes of size 1/A. Unfortunately, to leading order in 
A/kp, the mass shift went the wrong way. Presumably higher order terms (which are surely 
needed since the point flux limit is A/kp — > oo) will fix this. 

Once the kinetic energy is quenched all that is left of the i = j term is 

(141) 
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where we have emphasized that B* (a constant so far) is to be evaluated at rj as per our 
calculation. Suppose we change the external field by 5 A. Repeating the analysis with this 



field will lead back to the above formula, but with an extra piece per particle. This 



2m 



result also follows if SB varies very very slowly in space, say between one galaxy and the 
next. This is exactly the coupling of a magnetic moment ej 2m, mandated by SSH^. We did 
not have to attach it by hand, it emerged naturally. Note that the m in this term did not 
get renormalized. We could also write this term as follows 



? = I £x ^ p{x) v x6A = ~I £x6A ' j ™ 9 (142) 



e5B (r 

where 



g 

\mag = — Z X Vn (143) 

2m 

is the current associated with uncancelled cyclotronic currentd^E. It follows that if we add 
such a 5 A, it will lead to a charge density 5p = K 00 ^-- Taking the ratio of the applied 
vector potential to the induced Sp, we obtain 

K 01 = %q^-K m . (144) 
2m 



Note that we did not calculate K 00 , but merely showed that Kq\ has a piece proprtional to 
it with a factor ^ . 

Let us now turn to the i ^ j pieces in Eqn. (|135|) We treat the q sum as follows: 



^2 ex P -iqiri - rj) = 5 2 (ri - rj) - ^ ex P -^( r » - r j) ( 14 5) 

g q>Q 

The delta function vanishes on spinless fermion wavefunctions or on hard-core boson 
wavef unctions, and can safely be dropped. (Note that when H acts on a wavefunction to 
its right, we can move this delta function through the H+ to get the desired zero since the 
exchange produces corrections of higher order in q.) We have thus reduced the i ^ j pieces 
to a short-range interaction to be lumped with (the RPA part of) Hs a , the contribution from 
the short range CS field <5a that was not made dynamical. What is their collective role? We 
know that they cannot feed back to our zeroth order hamiltonian since a vanishing kinetic 
energy for LLL fermions is an exact result in the noninteracting theory. On the other hand, 
in the case of v = 1/2, we can show that these terms affect the large q sector in the following 
very desirable way. In RPA, the large q fields (including the term we shifted from small to 
large q ) conspire to produce the plasmon pole with the right position and residue in the 
density-density correlation, as shown in Appendix A. The issue is the following: There are 
no oscillators for q > Q. The fermion mass has been renormalized to m*, but they still have 
to produce the cyclotron pole at eB/m at large q via a collective mode. One therefore needs 
a Landau parameter f% to restore the correct cyclotron pole@@. The i ^ j interactions 
(when shifted to large q) fulfil exactly that role. It is very satisfying that although our 
oscillators no longer exist for q > Q, the physics (plasmon pole and residue to order q 2 ) is 
continuous across this border. 
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A. Electronic charge density 



Consider next Eqn.( pT38|) for the charge density 

p0ld = gcopAo + At + _ ig g x , 

V 87T 1 + COS yUA ~^ 

and the operator that is equivalent in the physical subspace 

s-ga . [(AW + , t( _ ?))] _ x n , e -., (147) 

In the MR, any combination of the form 

7/ d (g) + (l-7)^^ (148) 

is an acceptable definition of the physical charge density. (For example, HLR use this 
freedom to write the Coulomb interaction entirely in terms of a.) We could canonically 
transform any such combination to represent charge density in the restricted physical space. 
In an exact calculation there will be nothing to choose between them. In our approximation 
scheme, however, the following combination stands out: 

preferred _ \ old ( ) , 2 P gg^Jg) (ug) 

^ " 2p + l P {q> + 2p+l 4tt 1 } 

= -|=cosMo(A(g) + A\-q)) + - iig(J> x Uj)e^). (150) 



3 

To see what is special about it, consider p, its projection to the LLL, which clearly corre- 
sponds to the ground state of the oscillators, and is obtained by dropping the terms linear 
in A and A^: 

p(q) = ^ - «2(E(9 x ",)' ""' )• (i5i) 

Now we know from the work of Girvin and JacfS and GMP0 that the projected density 
must obey the algebra of magnetic translations: 

[p(<i),p(<i')} = ^o(<ix<i')p(q + <i')- (152) 

in the small q truncation of the structure constant, a truncation which satisfies the Jacobi 
identity )H. Our p obeys this algebra. There are two points that should be noted in this 
connection. 

• The magnetic algebra is expected upon projection to the LLL only in the original 
electronic problem, i.e., between physical states. In evaluating a product or commuta- 
tor, we must therefore restrict intermediate states to be physical. The exceptions are 
operators that do not mix physical and unphysical states, i.e., gauge invariant ones. 
These can be multiplied freely in the full space. Since we computed the p commutator 
with no restriction on intermediate states, we expect that p is gauge invariant^ and 
that we can freely multiply our p of Eqn (|150|) in the enlarged space without worrying 
about constraints. 



26 



• The algebra contains terms up to cubic order in q in the right hand side. These can 
arise from the terms we have kept and terms we have neglected. The truncation we 
have chosen has the virtue that the algebra is satisfied with the terms we have, with 
no need to appeal to higher order terms. Presumably the leading higher order terms 
cancel out in this preferred combination. 

For these reasons, we shall use this p(q) and the corresponding p old : 

P° ld = ^^(A(g) + At(-g)) + ^^e-^-^E(gxn,)e-^ (153) 

when we need a formula for the electronic charge density. Since the particle kinetic energy 
has vanished, the entire hamiltonian will be given by the Coulomb interaction written in 
terms of the above density. It will evidently be gauge invariant. 

Our formula for the charge tells us that the CF has charge 1/ (2p + 1). When p = oo i.e., 
v = 1/2 , this part of the charge vanishes and we are left with just the second dipolar piece. 
The latter couples to an external scalar potential $ as per 

H ext = (-e)Ep(-?)%) = -iel 2 J2J2(<lXPjy qXj ®(<l) = E- e ^x^-E(x,)(154) 

1 1 j j 

where E is the electric field. Note the anticipated value of the dipole momentSSS. 

We have argued that our charge density is gauge invariant in that it does not mix physical 
and unphysical states. [] This means it must either commute with the constraints or have a 
commutator proportional to them. Indeed we find 

[p(q),x(q r )}=u 2 o(qxq')x(q + q')- (155) 

Our joy is tempered by the algebra of constraints. We find 



[x(q),x(q ; )} = H 



2 (q' x g) 

COS 2 /J,\q 



P(q + q') + n - r r 2 £((<? + q') x U^e****' 
(1 + cos /iAo) 2 ~ 



.(156) 



Note that the operator in the right hand side is not quite the constraint; the term of order 
(q x q')(q + q') has the wrong coefficient. In contrast to the charge algebra which closed and 
was correct even to order (q x q') (q + q') , the constraint algebra closes only to leading order 
if we use the truncated expressions. Thus we have a situation where the Faddeev-Popov 
method applies ( the hamiltonian commutes with the constraints to give constraints which 
in turn form an algebra) only to leading order. 

What if we use the projected version of Eqn. (|146|) for the charge? The reader may verify 
that GMP algebra is satisfied and the constraints commute with the hamiltonian, both to 



Usually we say that an operator is gauge invariant if it commutes with the constraints, which 
are the generators of the gauge transformations. However if we merely require that the operator 
not turn physical to unphysical states, it is sufficient if its commutator with the constraints is 
proportional to the constraints. The Faddeev-PopovH method applies to this general case. 
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leading order (q x q'). In this version the charge of the CF will not be transparent, but the 
conservation laws and Ward identities will be. These comments apply to a models obtained 
by truncating operators to the orders calculated. Non-RPA and higher order terms in q can 
modify the commutators0. We do not pursue this issue for the gapped states for which we 
will propose an approximation scheme that ignores the constraints. 

Of all the results that flowed from the canonical transformation, only 1 / m* = requires 
Q = kp. For example the CF charge (l/(2p+ 1) ) and dipole moment are true for any 

B. Hall Conductance 

A common question one has for v = 1/2 is who carries the Hall current if the CF is 
neutral. The answer, for general v = p/(2p + 1), and not just for v = 1/2, is that the 
oscillators carry the entire Hall current in the clean limit considered here. We show this as 
follows. 

• We couple the system to an external potential <3>(g) using the formula Eqn.([L53|) for 

pOld 

• We recall that in the passage to the FR, the transport current due to particle motion 
got exactly canceled so that finally : 

Jt%) = ^ph A{q) . (157) 

V Z7T 

• In view of the above, we focus on just the oscillator sector, find (A) due to the potential 
and from it, (J). 

So we begin with 

H osc = ^2u c A\q)A(q) - e E <$>{q)p° l \-q) (158) 
q q 

= Y,"cA\q)A(q) - eJ2®(q)-7={A{-q) + A\q)) cos ^A (159) 



in 



E 



, + gecosuA(, / , 
A\q) - ^ ° $(-g) 

:7T0J r 



gecosuAo,/ x 



+const. (160) 
(161) 



The new ground state of the oscillators is found by shifting A as follows: 

(A(g)) = ^|^cos / iAo$(g) (162) 



The ground state electromagnetic current is 

cos /zAo 

7 2T < '™ ' h 



((-e)jf(q)) = (~e f +UJc < A(q) >= -^q + u*(q) (163) 
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where we have reinstated % = 1 = h/2ir to show that we have the correct value a xy = ve 2 /h. 
By translation invariance (in the clean system) the Hall conductance has to have this value. 
What is significant is that it all comes from the oscillators. In general we can write in our 
approximation scheme the following expression for longitudinal and Hall conductivities 

_ ^osc i ^particles f-\dA\ 

since the charge p old and hamiltonian separates into pieces, one for the oscillators and one 
for the particles, as shown in our letter^. (The particle part happened to be negligible 
for the Hall conductance in the clean system. When we turn on interactions, the particle 
current will acquire a piece proprotional to e 2 . ) This additivity of „, the has also been 
independently obtained by D.H. Leel, who statrted with the bosonic theory. In contrast, in 
the usual CF fermion theories in which the CF is fully charged, one adds resistivities. The 
concerns of Lee et a@ on the value of therefore do not therefore apply to our approach. 



VI. INTRODUCING INTERACTIONS 

So far we have been considering just the noninteracting theory. A unique ground state 
appropriate to an interacting theory emerged in the MR upon keeping just Hq and dropping 
Hj and H n . This led to good wavef unctions, but bad dispersion relations for the oscillators 
(wrong frequency at q = except when A* = 0) and a particle kinetic energy of order 
1/m instead of zero. We then dealt with Hj and Hjj: eliminating Hj by the canonical 
transformation and dropping Hjj as non-RPA. This led to oscillators of the right frequency 
and particles with no kinetic energy and just a magnetic moment e/2m which coupled to 
variations in the background B field. These are the correct results for the noninteracting 
case. 

We are now going to include interactions, illustrating the procedure with a Coulomb 
interaction cut off in the ultraviolet at q = Q so that it may be described entirely by our 
oscillators. (This rounding off of the Coulomb potential at short distances may be a realistic 
description of the experimental system with finite thickness.) The full hamiltonian is: 

H T = E"cA\ q )A(q)+j: e ^ + lf: P ° l %- q ) 2 ^p°%) (165) 

q i Z?T1 Z q V 

p ° ld{q) = 9 ^W L{Aiq) + AH ~ q)) + ^TT - u ^ iq x n >~^)- ( 166 ) 

When we expand out p old (—q) v(q) p old (q) we will get a piece that renormalizes the os- 
cillator frequency to order e 2 q (which for short range interactions would have been order 
e 2 q 2 v(q)) as per Kallin and Halperin@), a mixing of order e 2 q between oscillators and par- 
ticles and a p{—q) ]1 ^-p{.q) term in the particle sector. The mixing can be eliminated by a 
further canonical transformation, which can be done to lowest order in e 2 q. It will modify 
all previous formulas for p old , J and so on to terms of higher order. We ignore these changes 
in our analysis, which is limited to lowest nontrivial order in q and e 2 . 
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Before going to the particle sector for a detailed analysis let us note that if in the density- 
density correlation K pp we consider just the oscillator part of p old : 

p M ^ q -^(A( q ) + A\- q )) (167) 



■57T 

we find the cyclotron pole 

v , x q 2 uj c cos 2 p\ . . 

47r(u^ - (lj + irjY) 

whose residue obeys the sum rule 

oo q^rni 

ImK 00 (uj)ujdu) = . (169) 

o 2m 

Since we work in unit volume n = N, the particle number. 
Turning to the particle sector, the hamiltonian is just 

Q ne 2 

H part = J2p(-q)—p(q). (170) 

q Q 

(We have dropped the eB* /2m term assuming that the external field is homogeneous). 

Recall that GMP§ be gan with the above hamiltonian, the LLL assumption, and the 
magnetic algebra of the p's, and proceeded to derive the spectrum of collective excitations 
in the Single Mode Approximation. What is new here, as far as low energy physics is 
concerned? The first point is that we have not just the algebra, but a concrete realization 
of the algebra in terms of canonical variables in terms of which we can try to do explicit 
calculations. This is similar in spirit to the recent work of Haldane and PasquieiS who 
started with an algebraic formulation of the LLL problem and then switched to a concrete 
realization of the algebra in terms of canonical operators. Next, our approach derives the 
representation from first principles, and shows how a field theory in the full Hilbert space 
leads to LLL physics upon identifying the plasma oscillators and then freezing them. Finally 
there are many low energy responses depending on m, which we obtain by starting in the 
full space and systematically separating high and low energy degrees of freedom. 

The particle hamiltonian, written out fully, (dropping the e5B/2m term) is 



Hpart 




1 i,j Q 



Q p 2 iqixi-Xj) -72 Q _p2 

vr^e + _^o_y^Y:{ — (qxU j )e^-^-h.c.} (172) 



q 1,3 t V r 1 -V "f i q 



This hamiltonian affords a concrete, microscopic, and detailed realization of composite 
fermions. Not only does it describe particles that see the weaker field B* that is just right 
to fill p Landau levels, it also describes particles with local charge e* = e/(2p +1). The 
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main attraction of Jain's CF picture is that one can think of the original problem in terms 
of roughly noninteracting particles. Since our description shows the right charge and a mass 
of Coulombic origin at tree level, we may expect that only small corrections separate what 
we see in the hamiltonian from what finally happens. 

The hamiltonian is, however, not that of free fermions in the weaker field. Indeed, this is 
not to be expected, because the Hilbert space of free fermions is too big to describe the LLL 
physics that is left after we have frozen the oscillators. That is why we have constraints to 
reduce the size of the Hilbert space and the hamiltonian has to be gauge invariant under the 
transformations they generate (to the approximation we are working in). The free particle 
hamiltonian does not have this property. 

However H part contains within it a free particle piece which could serve as a starting 
point for an approximation. It is given by the i = j part of the dipole-dipole term 



Tie 2 



H£% = (g 2 |n/sinXn) (173) 

q j Q 

in -12 

^ 2m* v ; 

3 

m* q q 6 

where we have used the fact that Q 2 = Ann = eB cos 2 /i\ and cos 2 u\ = 2p/ {2p + 1) = 2v. 
Notice how the interaction has generated an effective 1/m*. Its exact value is not universal 
and depends on the rounded Coulomb interaction employed. 

Read has argued0 that the origin of the kinetic energy is the binding of the CF to the 
correlation hole, and that the effective mass is the curvature of this potential energyreadl. 
(The curvature in real space gives the mass because the momentum of the fermion is a 
measure of its dipole moment.) If we consider the potential we have used, (the Coulomb 
interaction cut off at Q in momentum space, so that it is rounded off to a finite value at 
r = 0), and calculate its curvature at r = 0, we find that it is given precisely by the above 
formula. Of course the bare potential gives just the bare mass which can be renormalized bv 
the other terms in H part . It is encouraging that numerical work of Morf and d'Ambrumenilej 
gives 1/m* = .2c 2 Iq for the untruncated Coulomb interaction at v — 1/2 compared to our 
e 2 Zo/6 for the truncated Coulomb. In any event, the aim of the above calculation is merely 
to show that we have a means of computing 1/m* in a zeroth order approximation and that 
it has the right size and physical origin. To our knowledge our earlier wor was the first 
instance this was done. The exact value for a given potential is nonuniversal and dependent 
on small and large g's and beyond the reach of our theory. 

Note that the kinetic energy that comes from the i = j term does not constitute the elec- 
tron's self-energy, but rather the interaction energy between the electron and its correlation 
hole produced by the (deficit of) other electrons. One may ask why we did not banish the 
i = j term of the Coulomb interaction in the MR from the outset. The problem is that 
we only know how to canonically transform the full charge density, summed over all parti- 
cles. This density obeys the magnetic algebra and functions of it are gauge invariant. So 
our strategy has been to add the i — j term in the MR, where it constitutes an innocuous 
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chemical potential shift, work with the product of the full densities in constructing the in- 
teraction hamiltonian and then transform the latter. The magnitude of the kinetic energy 
so generated fits our expectations based on Read's arguments and points to the soundness 
of the procedure. 

Returning to general fractions, we have to understand why the CF are weakly interacting 
and may be described approximately by H p ^ e , for this is essential to the success of the 
CF picture of Jain!. The smallness of e 2 does not matter since it is a prefactor to the 
whole hamiltonian. Looking at the neglected terms, we see two possible small dimensionless 
parameters: 1/p and ql . Our preliminary analysis of matrix elements (relying on the work 
of work of Dai et and Chen et a&) shows that there is huge window (which gets larger as 
v — > 1/2) where the dipole can dominate the monopole because ql^^/p * s l ar g e even though 
qlo is small, Zq being the magnetic length appropriate to the reduced field A*. Not having 
completed the detailed calculation we merely outline a possible strategy. Suppose we want 
Koo(q), the retarded p old — p old correlation function. Since p old does not link physical to 
unphysical states, its two point function in a gauge invariant ground state can be calculated 
without regard to the constraint on intermediate states. However we do not have access to 
this gauge invariant ground state. We must argue that the numerical value of the correlation 
will not be too different in an approximation to the ground state. For the latter we may 
use the p- filled Landau level ground state of H^rt- We could also include the neglected 
interactions terms in Eqn. (|172|) in an RPA calculation. We are in the process of evaluating 
this strategy by computing some correlations that are known numerically or experimentally. 
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A. v = 1/2 

We now focus on v = 1/2 to which a great deal of attention has been given following the 
work of HLR. They took the notion of a Fermi surface at v = 1/2 seriously and deduced 
many experimental consequences which have been verified^. 

HLR calculated numerous response functions in the RPA working in the CS representa- 
tion we began with (before the introduction of the oscillator coordinates). Let us focus on 
the following results£j: 

oo 

ImK 00 co dco — q^ (176) 

oo 

ImK 00 duj ~ g 3 logg (177) 

o 

roo 

/ ImK 00 co" 1 dco ~ 1/v (q) = q (for Coulomb) (178) 
Jo 

where the plasmon contribution has been removed. 

What does our formalismpredict for this response function? The effective hamiltonian, 
the constraint x(l) an d p arella: 

„2 ;4 Q o 2 

H P art = E ^ + j E E ^(?x ft )( 9 x ft )e^-^ (179) 
X(q) = E e~ im +yE(? x Pi)^ = (180) 
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l 2 J2(q x Pi) (181) 



To calculate K 00 , we need an approximation scheme. One possibility we explored starts 
with the p 2 /2m* term in Eqn.( |179|) and includes the other terms and constraint in the RPA. 



While this appears reasonable and familiar, there are some dangers in this CcLSG, clS will be 
explained shortly. The result in the small q limit is , 



-ft 00 



4tt{uj 2 - {uj + if)) 2 ) 



771*1 o „, o , /—r t ox ^ UJ 



2tt l 2 v 7 J gw* 

= -ft"oo,osc + Koo,cf (182) 

We find that the answer is not changed (in the small q limit) as we change the constraint 

from Ei e~ iqn + f Y.M x Pi)e"^ = to Ei e"^' 1 = 0. 

The above formula includes the plasmon contribution as well. Note that the CF contri- 
bution has no structure other than the particle- hole branch cut. It follows from the above 
that 
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ImK CF , m u a dw = fcrVTV** - ^,0^ io\ « = 0, ±1 (183) 
o H-s/7TL {b/2 + a/2) 

The ratio of the a = 1 and a = integrals, which gives the average frequency of the 
excitations, goes as 

< u >~ q. (184) 

The q dependence of the first two sum rules agrees with that of HLR (up to logarithms) but 
the last is quite different: in our case the compressibility vanishes as q 2 and is independent 
of the interaction, while theirs goes as l/v(q) ~ q for the Coulomb case. In the last sum rule 
very small cj's dominate. In HLR, this region is controlled by the overdamped mode, which 
goes as uj ~ iq 3 v(q). This mode arises due to magic gauge cancellations in the infrared. 

We did not see such a mode or such magic cancellations. This could be due to an 
improper treatment of the constraints by our RPA. Could it be that if the constraints are 
treated properly, negative powers of q appear in the correlation function and neutralize the 
q 2 coming from the dipolar nature of charge? 

Let us look at our constraints and the gauge transform they generate. While we do not 
know them exactly, there is no doubt that as q — » 0, they reduce to 

x(q) = J2^ qri=0 - (185) 

Their action on the particle coordinate and momenta are as follows: 

Xj -> rj Pj -> Pj + qe- iqr i (186) 

Thus x essentially shifts the momenta of the particles. This means that the Fermi circle 
cannot be a ground state- it can be moved around by a gauge transformation. This is 
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reminiscent of Haldane 's finding that the particles like to cluster around each other in mo- 
mentum space, but with no preferred origin. In our analysis this reflects a subset of the 
gauge symmetries of the effective hamiltonian. 

Now, in any gauge theory one resorts to gauge fixing by choosing an operator whose 
commutator with x is non zero and freezing its value, say to zero. For example, in QED, the 
constraint V ■ E = is accompanied by a gauge fixing condition V ■ A = 0. Such gauge fixing 
kills all spurious fluctuations. In our case it may justify our fixing the Fermi sea arbitrarily 
at the origin. On the other hand, in a gauge theory, even after gauge fixing, there can be 
genuine gauge invariant fluctuations of arbitrarily low energy. We did not see them in our 
treatment, and the question is whether they would arise in a different approximation and 
if so, produce inverse powers of q that lead to a finite static compressibility for short range 
forces and and one that vanishes as q for the Coulomb case. 

Recently Halperin and Sterni have given an existence proof that such a thing can hap- 
pen, and indeed happens in a model that arises within our formalism. Consider the nonin- 
teracting hamiltonian Eqn.( |135|) , drop the oscillator part and the magnetic moment part, 
and set v = 1/2 so that U = p. Then 

^£f|-^££I>-^-^ (187, 

i % j q 

In the Q = limit, the hamiltonian is invariant under a uniform shift of all momenta. For 
this reason choose Q very small. The formula for charge is still dipolar since that does not 
depend on Q. Halperin and Stern compute the density- density correlation in this model and 
show that inverse powers of q indeed appear due to gauge cancellations and modify render 
the system compressible. 

The authors do not claim that this is the problem we set out to solve - it has explicit 1/m 
dependence, no Coulomb interactions, describes particles that are bound to fat flux tubes 
and do not obey any particular statistics except outside a distance 1/Q. The role of this 
calculation, as stated above, is to demonstrate that inverse powers of q can and do appear 
in a model with the same gauge symmetries, at least as q — > and that dipolar particles can 
be compressible. The authors however state that the same conclusion applies to the model 
at hand and promise details soon. This result was also announced by D.H. Lee in a revised 
version of Ref . (|4*2"| ) . 

Note that even if we regain the HLR formulas for correlations, the physical picture is 
now much more satisfactory. We have dipolar fermions and a special hamiltonian whose 
symmetries allow for compressibility despite the dipolar coupling. This hamiltonian and 
its symmetries are not postulated but derived from the microscopic theory. We also have 
oscillators which contribute additively to Hall conductance to complete the picture. Thus 
we reconcile many seemingly contradictory properties of the CF's. 

One can seek other approximations besides our RPA and see if such cancellations occur. 
Another possibility is to use Eqn. ( |146| ) for charge so that 



3 1 j 



(188) 
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In this case the constraint 



X 



3 3 



(189) 



commutes with the charge (and hence the hamiltonian) to leading order and gauge invariance 
will be easier to implement. 

In the meantime we can ask if there are kinematical regions where our RPA analysis 
may hold independent of how well it treats the constraints. The first two moments of Kqq 
and intuitive arguments suggest that higher the frequency, the better the prospects. With 
this in mind we compared our results to the surface acoustic waves (SAW) experiments of 
Willett et al (1993) in the frequency range of ~ 3.7 — 6.1GHz. 

The fractional shift in the velocity of SAW, as explained to us by SimonJ^ is given by 

Sv./v, = 3.2 ■ 1(T 4 [1 - to(v(q) K 0Q (q, u))] (190) 

where v(q) = 2ire 2 /q and 3.2 ■ 10~ 4 is a material dependent constant referred to as a 2 /2 in 
the literature and 3? denotes the real part. The formula is set up so that Sv s = for a free 
fermion system in zero field, for which dtKoo(q,uj) = l/v(q). 
Using our result for K Q0 we obtain 

c )ttp 2 m* 

5v s /v s = 3.2 • 1(T 4 (1 ^{qhf). (191) 

eq An 

Feeding in 

_L_ Ce 2 l 
m* e 

and 



(192) 



u> 2tt ■ 10 9 / 

q = — = 193 

V s V s 

where / is the frequency in GHz, we get 

10 4 • 5v s /v s = 3.2 - (194) 

upon using v s = 3000m/s, and k,F = 1/Iq = 93/im -1 . Notice that the dielectric constant 
cancelled out. 

Extracting from Fig. 29, of Willett's 1997 review, the points 

(/, 10 4 ■ Sv s /v s ) = (3.7, 1.85), (4.3, 1.66), (5.4, 1.37), (6.1, 1.17) 

we obtained a least square fit 

(Svs/vs) ■ 10 4 ~ 2.9 - .28/. (195) 

Upon comparing to Eqn.( |l94] ), the slope translates to a value of C ~ 0.13. (For reference, 
the zeroth order value for the rounded Coulomb potential was C ~ 0.17.) The intercept, 
given by the on the material dependent constant 3.2 • 10~ 4 is within ~ 10% of the data. 

We are aware that comparing theory to experiment (rather than to another theory) is 
lot more involved. We do not claim to understand in any depth why our RPA should work 
in this kinematical region. But we are intrigued by this agreement with experiment and 
present it as such. We propose to study other response functions in this kinematical region 
to pursue this question. 
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VII. CONCLUSIONS 



The aim of this paper was to start with the problem of planar electrons of mass m, charge 
— e in a perpendicular magnetic field —B and see how far we could go towards understanding 
the FQHE within some approximation scheme. We discussed the fractions v = p/(2ps + 1) 
focusing on s = 1, the extension to other values being very direct. 

We began by following the standard procedure of attaching point flux tubes to go from 
the electronic to the CS representation. In the hamiltonian description this led to the 
introduction of the nondynamical CS field (Vx) _1 27r/p when I quanta were attached. The 
average value of this (due to the average charge) cancelled some or all of the external field A 
leaving behind an A* which either vanished for v — 1/2 or v — l/(2m+l); or for v = p/2p+l 
was just right to fill p Landau level a la Jain. 

At this point we took the crucial step of enlarging the Hilbert space to include the 
canonical pair a(q), P(q) which described transverse and longitudinal vector potentials with 
< q < Q = kp. Thus the number of additional degrees of freedom equaled the number 
of electrons n. This led to n constraints on physical states (a(q)\physical >= 0). We 
made a further canonical transformation that got rid of the nasty field (Vx) _1 27r/p and 
changed the constraint to qa = 2nlp. The hamiltonian in this middle representation (MR) 
did not contain any dependent fields. We found that if we dropped the terms Hj and Hjj we 
obtained a solvable model of n plasma oscillators described by a, P and n particles, either 
free or in a field just right to fill p Landau Levels. The ground state (product) wavefunction 
for this case, duly projected to the physical sector using the constraint, gave us the well 
known correlated wavefunctions of Jain and Rezayi-Read, (except for the overall projection 
to the LLL). The oscillator part of the product was the Jastrow factor associated with flux 
attachment, but it came with the magnitude and phase of the correlation zeros, as well as 
the gaussian factors with right magnetic lengths. 

We identified operators that created Laughlin's quasiholes. The corresponding wavefunc- 
tions had the right Gaussian factors for the quasiholes, and the normalization factor from 
which one could infer their statistics!. We constructed the operator that could be identi- 
fied with the Read order parameter, but pointed out some limitations of this identification. 
Our composite boson and fermion operators implemented the physical picture extracted by 
studying Laughlin's wavefunction- that electrons like to bind to zeroes of the wave function. 
Only the zeros produce a charge deficit to which the electron is attracted; it is not attracted to 
just flux, i.e., the phase of the zeroes. It appeared as if every nondegenerate ground state of 
the interacting theory had a counterpart in the free theory of oscillators and particles upon 
projection. All these results were obtained before turning on any explicit e — e interaction 
although it was understood that such an interactions exists in order to stabilize the unique 
ground state. 

There were however some problems with the MR that did not surface in the quest 
for wavefunctions. The plasmons had the wrong frequency (unless A* vanished) and the 
particles had kinetic energy of order 1/m, when they were supposed to have none (in the 
noninteracting case) or one of order 1/m* ~ e 2 l$ in the presence of interactions. To cure 
this, we performed yet another transformation to the final representation (FR) to decouple 
the oscillators and particles. This was done in the small q limit and within an RPA applied 
to the final variables!!. The result was as follows. 
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First, the constraints were entirely in the particle sector. Denying n particles in two 
dimensions n degrees of freedom was just right to project out the LLL physics. The particles 
and oscillators separated in the hamiltonian. The particle energy completely vanished except 
for the coupling of a magnetic moment e/2m to the external field, a coupling whose existence 
had been predicted by SSH. Thus the desired quenching of particle KE, i.e., 1/m* — was 
accomplished. The oscillators ended up with the correct frequency. The formula for p old , 
the charge density of the electrons, contained three pieces: one linear in the oscillators, 
one proportional to the particle density with a charge l/(2p + 1) and a dipole piece. Thus 
we regained the correct charge of the CF and the dipolar nature of the charge for v = 
1/2 as anticipated by Read. If we dropped the oscillator part of the charge (valid in the 
oscillator ground state sector) we obtained a projected charge p which obeyed the magnetic 
translation algebra of LLL charge density first noted by GMP@ (in the small q limit). 
This projected density was gauge invariant i.e., did not mix physical and unphysical states, 
and could be multiplied freely in the enlarged space. Charge-charge correlations and Hall 
conductance were additive over the oscillators and particles since the formulae for charge 
and hamiltonian were. The oscillators were found carry the full Hall current in the clean 
limit we are considering. 

We did not try to obtain improved wavefunctions from the FR since the transformations 
back to the electronic wavefunction were no longer simple. This simply corresponds to the 
fact that projection to the LLL (accomplished by the passage to the FR in our scheme and 
by the V operator by Jain and Rezayi-Read) does something very complicated to the wave- 
functions we derived in Section IV. Since we went from the MR to FR in first quantization, 
our main results- quenching of kinetic energy and the reappearance of an interaction de- 
pendent 1/m*, the value of e*, the constraints, additivity of response function over particles 
and oscillators- all apply to the bosonic case (u = l/(2s + 1)) as well. Also all the results 
in the FR, except for the cancellation of 1/m in the noninteracting case, were insensitive to 
the choice Q = 

We turned on the Coulomb interaction and found that the particle sector now contained 
a kinetic energy term. A zeroth order value for 1/m* was easily found. It had the right 
order of magnitude and physical origin. Our hamiltonian had a gauge symmetry generated 
by the constraints. In the infrared limit, these transformations were seen to be a translation 
in momentum space, reminiscent of the drifting Fermi sea seen by Haldane. We proposed a 
computational scheme that was quite reliable for gapped states. For v = 1/2 we described 
an RPA calculation of K 00 , the retarded 

p old _ p old 

correlator. Two truncations of the con- 
straints gave the same answer. Two of the Kqq moments were the same as that of HLR 
(up to logarithms) while the compressibility vanished as q 2 in our case versus approaching 
l/v(q). We pointed out that the compressibility formula could be altered by gauge cancella- 
tions that we did not find, cancellations that have been shown to occur in a related model by 
Halperin and Sternlll. We emphasized that even if in the one ends up with the HLR correla- 
tion functions, our method has facilitated the reconciliation many seemingly contradictory 
properties of the composite fermion such as its dipolar charge and compressibility. 

We went on and compared our results to the SAW experiments of Willett et al on the 
expectation that perhaps in that kinematical region the precise treatment of constraints 
would be unimportant. We presented the agreement with data as an intriguing fact that 
needs to be understood. 
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In summary, we have proposed a way to extract a low energy theory for the FQH states. 
Composite bosons and fermions with the frequently quoted properties arise at tree level 
in the FR, but one finds that the oscillators are needed for a complete description. The 
most important ingrdients are explicit expressions for the hamiltonian and charge-current 
operators of the CF and CB. Our description relies on a small q and RPA approximation 
(in the final variables), the precise nature of which we do not fully understand. The same 
applies to the gauge constraints or gauge symmetry in the FR, of which we are sure only of 
the extreme small q limit. As we, along with others , continue to work on clarifying these 
issues, we have provided here many details of our machinery in the hope that they may be 
fruitfully employed by a wider audience. 
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APPENDIX A 



We show here that the short distance terms H$ a along with the i ^ j terms in the third 
sum in Eqn.(|135]) that we shifted to large q with a change of sign, conspire to produce the 



plasmon pole of correct location and residue in the case v = 1/2, when II = p. Let us focus 
on just this part of the hamiltonian, again keeping only zeroth order terms in q. (Note that 
we do not claim q in this region is small, only that we are working consistently to a given 
order in q.) Recalling that in this region p = p, 

1 
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EEE^'Pj e 



-iq(n-rj) 



+ E E |e^a(g) + £ *a(g)*a(-g) (196) 
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We now switch to the functional formalism. We perform a Hubbard-Stratonovich transfor- 
mation of the first term using a field V = iqVi — iz x qVr- The lagrangian density is (upon 
dropping the 5 symbols in front of the fields, all of which have q > Q) 
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~a(q)a(-q) + ^(V* + V 2 L ) + a (-q)j-a(q) 



(197) 



It is now straightforward to calculate the inverse RPA matrix propagator for the 4- 
component field D = (a , a, Vl, Vt) = (-Do, D t , Dl, Dt) using 

D- 1 = + U° (198) 
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where II is the free-field response function. In particular 



[n°] 00 (g^) 



2n d6 k F q cos 8 
47r 2 u) — v*q cos 9 + irj 



(199) 



while the TT (or tt) , LL, and OL/LO components differ by the additional factors of 
(hp/m) 2 sin 2 9, (/cf/tb) 2 cos 2 9 and ±i(kp/m) cos9 respectively in the numerator of the in- 
tegrand. We do not list these integrals; they may be found (for the time-ordered instead 
of retarded functions) in the Appendix of Reference (0). Notice that even though 1/m* 
(or v*) vanishes here, we keep it in the denominator since the answers can have m* in the 
numerator. Also since we are working to leading order in q, we do not keep the q 2 /2m* in 
the energy denominator. These functions are needed only in the limit 
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The inverse matrix in this limit is 
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(201) 



It is now straightforward to verify that as x —>■ oo (which is the only region we have here) 
this propagator has only a plasmon pole at c^o- The density- density correlation, related to 
Dhy ' 



K = K — K DK 



(202) 



(where K° is the noninteracting limit) has the following matrix element: 



00 



47T L) 2 



{uj + ief 



which obeys the f-sum rule 



q 2 nn 

ImK(uj)ujduj = . 

o V ; 2m 



(203) 



(204) 



It is important to note that if we had not shifted the i ^ j pieces from small to large 
q the plasma pole would have ended up at Ann / V 'mm* . This term would have also been 
unwelcome at small q since the fermions are to have zero hamiltonian. All in all, the shifting 
from small to large q seems decidedly the right thing to do. 
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